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/. Introduction 

WHILE making a comparative study of the scattering of light in optical glasses 
with the incident light in different states of polarisation, Krishnan (1936) 
found that the depolarisation factors p, and p, corresponding to incident 
unpolarised light and light polarised with vibrations vertical respectively, were 
of the same order of magnitude as those usually observed in ordinary gases. 
On the other hand, when the incident light was polarised with vibrations 
horizontal, the depolarisation factor p, instead of being equal to | as in the 
case of a gas or a liquid, was less than | as in the case of colloids containing 
particles not small compared with the wave-length of light. This anomalous 
depolarisation was explained by Krishnan as due to the existence of mole- 
cular aggregates in glasses. The actual magnitude of the effect is rather 
small and hence it is difficult to observe visually whereas it is easily seen 
with many ordinary colloids. The observation of the effect is also rendered 
more difficult by reason of the feeble intensity of scattering when the inci- 
dent beam is polarised with the vibrations horizontal and by the weak 
fluorescence exhibited by the glasses. Krishnan (1938) therefore, devised an 
objective method based on photographic photometry for demonstrating 
the reversal of polarisation in optical glasses. In this method the transversely 
scattered light is photographed through a double-image prism using incident 
polarised light. Photographs of the two tracks were taken on the same 
negative for different orientations of the Nicol (polarising the incident beam) 
about the mean horizontal position. Microphotometric records of the 
negatives published in his paper show clearly the reversal of polarisation 
when the incident beam is polarised with vibrations horizontal. 

It is well known that the light scattered by any medium whose elements 
of heterogeneity are not small compared with the wave-length of light, will 
be elliptically polarised when the incident beam is linearly polarised in a 
direction oblique to the plane of scattering. Elliptic polarisation has been 
observed by Hariharan (1942) in the light scattered by liquid mixtures very 
near the critical solution temperature, in which molecular clusters of appre- 
ciable size exist. The detection of elliptic polarisation in the light scattered 
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by optical glasses under similar conditions would give us an entirely inde- 
pendent and conclusive evidence for the presence of molecular aggregates 
in such media. It is the object of this paper to present such evidence. 


2. Details of the Experiment 


The new technique for the study of light-scattering developed by 
Sir C. V. Raman (1941) has been used with success in the present investigation. 
This new method is based upon the use of a Nicol for polarising the 
incident beam in any desired azimuth and of a Babinet Compensator for 
observing the transversely scattered light. It enables one to demonstrate the 
reversal of polarisation or anomalous depolarisation and also to detect the 
existence of an elliptic polarisation in the scattered light. 


The optical glasses used for the present study were the same specimens 
with which Krishnan (1936) carried out his investigations. These glasses 
were absolutely free from strains and did not show any measurable 
birefringence. Of the seventeen glasses only nine were selected for the present 
investigation. Their serial numbers were 1, 3, 4, 5, 6, 7, 8. 9 and 16. 


Sunlight was reflected into a dark cabin by means of a heliostat and 
focussed by a Dallmeyer photographic lens of variable aperture. A 
rectangular glass cell containing a concentrated solution of alum was kept 
close to the lens in order to cut off the heat rays. A Nicol mounted on a 
graduated disc was placed in the path of the incident light before it entered 
the optical glass under investigation which was placed inside a wooden box. 
By rotating the Nicol the incident beam could be polarised with the direc- 
tion of vibration inclined at any desired angle with the vertical. The 
specimen of optical glass to be examined was suitably blackened on the out- 
side excepting for three windows, one for the entrance of light, one for its exit 
and the third for the observation of the transversely scattered light. To 
view the light transversely scattered by the specimen, a Babinet compensator 
was placed with its principal plane inclined at 45° to the plane of observa- 
tion. The scattered light then passed through another Nicol oriented in 
such a manner as to transmit vibrations which had their electric vector hori- 
zontal. An eye-piece was used for making visual observations. The eye- 
piece was removed and substituted by a camera for taking photographs of 
the pattern of fringes in the field of vision. 


3. Sources of Error 


The main source of error in the present experiment is the presence of 
inclusions or specks inside the glass, which give rise to an intense polarised 
scattering. A preliminary survey was made of the seventeen optical glasses 
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available in this Laboratory. A strong beam of sunlight was focussed into 
the glass prism and the track inside was examined visually along the direc- 
tion of the incident beam. Any speck or inclusion of colloidal dimensions 
would give rise to an intense and localised forward scattering and would 
thus make itself visible. Only those glasses were selected for the present 
investigation inside which the scattered track was of uniform intensity and 
colour throughout the entire section of the specimen. 


Another source of error is the weak fluorescence exhibited by the 
glasses. As was pointed out by Krishnan (1936), in order to eliminate the 
effect of fluorescence completely it is necessary to use an orange filter in the 
path of the incident beam. With the orange filter the intensity of scattering 
would be considerably reduced with the result that it might not be possible to 
photograph the fringe system within a reasonable time, especially when the 
incident light is polarised with vibrations horizontal. Therefore, in the present 
investigation a blue filter was used in the path of the scattered light. The 
width of the transmission band of the filter used was about 500 A.U. With 
this filter it was found that the vertical and horizontal component of the 
scattered light were almost of the same colour in all the glasses. 


A third source of error is the parasitic illumination. This was avoided 
by eliminating all extraneous light other than that passing through the glass. 
A series of diaphragms and screens in the path of the incident beam cut off 
all parasitic light and there were similar diaphragms on the observation side. 


To avoid errors arising from the finite angle of convergence of the 
incident beam, the aperture of the condensing lens was kept very small. 


4. Results 


The following general features were observed in all the glasses when 
the plane of polarisation of the incident beam was rotated. When the inci- 
dent light was polarised with its vibrations vertical. the fringes 
in the Babinet compensator appeared very brightly in the field of 
vision. The compensator used was of the fixed type and was provided with a 
vertical cross-wire in the field of vision. The cross-wire appeared on the central 
dark band, indicating thereby that the scattered light which was partially 
polarised under such conditions, contained an excess of vertical component. 
As the Nicol polarising the incident beam was slowly rotated, the visibility 
of the fringes as well as the intensity of the pattern decreased progres- 
sively. At a particular inclination of the plane of polarisation of the inci- 
dent beam, the fringes ceased to be visible and the field of vision was 
uniformly illuminated. For this setting of the Nicol in the path of the inci- 
dent beam, the vertical and horizontal components of the scattered light 
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were of equal intensity, i.e., the depolarisation factor pg was equal to |, 
# is the inclination of the plane of polarisation of the incident beam to the 
vertical. The value of 6 for which pg is equal to | can be calculated, using 
Krishnan’s formula* (1939), 


Po = (1+ tan?@’p,) | (tan?@+- 1 /p,). (1) 


The observed value of @ for which the fringes disappeared, was in reasonable 
agreement with the calculated value. 


Long before the fringe system disappeared, it was found that the pattern 
had shifted laterally and the cross-wire instead of appearing at the centre of 
the central dark band, appeared to have shifted towards the next bright band. 
The shift increased gradually as the polarising Nicol was rotated. Finally when 
the incident beam was polarised with vibrations horizontal. the cross-wire 
appeared at the centre of a bright band, showing thereby that the scattered 
light under such conditions was partially polarised with an excess of hori- 
zontal component, i.e., p, was less than |. The existence of the anomalous 
depolarisation in glasses is thus indisputably established. For those 
settings of the polarising Nicol for which the cross-wire appeared neither 
at the centre of a bright nor at the centre of a dark band, but at some inter- 
mediate position, the transversely scattered light was a mixture of unpolarised 
and elliptically polarised light. On rotating the polarising Nicol further 
through 90°, the features described above appeared in the reverse order. 
Typical photographs showing the general features of the phenomenon are 
reproduced in the accompanying plate. 


In order to measure the phase difference of the component vibrations 
in the elliptically polarised scattered light, the following procedure was 
adopted. The eye-piece used for visual observation was removed and the 
scattered light as seen through the compensator and the analysing Nicol was 
photographed. From the photograph the position of the cross-wire in 
relation to the nearest dark band was noted. Assuming the phase difference 
to be equal to = corresponding to a shift of the cross-wire from the dark 
fringe to the next bright fringe, the phase difference between the two compo- 
nents of the scattered light corresponding to this particular position of the 
cross-wire was evaluated. In this way the values of the phase difference 
for the various glasses for a series of orientations of the polarising Nicol 
were determined. The results obtained with glasses Nos. | and 5 are gra- 
phically represented in Fig. 1. 


* Perrin (1942) has since shown that Krishnan’s formula must be true for any symmetrical 
medium. 
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Table I gives the values of p,, p, and p, for the glasses examined. 
TABLE I 





Observed 


4 
2 | 





| ellipticity 
| is found 
| 


XP 





| 


26° 


| | 9° 
| *3 33° 


| 8° 
30° 
27° 
| 12° 
| | i 
| . 24° 

| | . 15° 

| } { 

The values of the depolarisation factors given in the table are those reported 
by one of the authors in an earlier paper (Krishnan, 1936). The column 6 in the 
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table gives the values of Ao, which is the difference between the observed 
values of p, and the anisotropic part of 9, which can be considered to be 
equal to 2o,/(1+ 9,) to a first approximation. Ap, therefore represents the 
depolarisation due to the finite size of the scattering centres. Ap, can be 
expressed in terms of either p, and p, or p, and p,, and the corresponding 
expressions are given below. Equation (2) or (3) can be used for the 


APpy = Py — 2Pp/(1 + py) (2) 


an =(- HC!) ° 


evaluation of Ap, depending upon the accuracy in the determination of p, 
or p, respectively. The last column in Table I gives the range of 4 over 
which the scattered light is elliptically polarised. 


5. Discussion 


From an analysis of the results obtained, ‘t is clear that in all the glasses 
examined, the light scattered in the transverse horizontal direction when 
the incident light is polarised with vibrations horizontal, is partially polarised 
with the horizontal component stronger than the vertical component. These 
results confirm Krishnan’s conclusion that the optical heterogeneity in glass 
cannot be identified as due to the individual molecules but is due to 
molecular aggregates, thus presenting an analogy with the case of liquid 
mixtures in the neighbourhood of the critical solution temperature and in 
protein solutions. 

According to the theory of Mie, in directions other than parallel and 
perpendicular to the incident electric vector, the light scattered by large 
spherical particles would be elliptically polarised. The existence of such 
elliptic polarisation in the scattered light has been detected experimentally 
in the case of emulsions by Darbara Singh (1942), in the case of liquid 
mixtures very near the critical solution temperature and in silver sols by 
Hariharan (1942) and in some typical colloids by Rao and Krishnan 
(1944). The present investigation reveals that in glasses also the scattered 
light is elliptically polarised when the incident light is plane polarised, 
The curves of ellipticity show a striking similarity with those obtained with 
colloids in which the particles are known to be not small compared with the 
wave-length of light. These observations also support the view that mole- 
cular aggregates do exist in amorphous substances like glass. 

Another important fact which emerges out of the present investigation 
is that those glasses which give larger values of Ap, exhibit ellipticity to a 
correspondingly larger extent. Compare columns 6 and 7 in Table I. This 
js not surprising because a larger value of Ap, corresponds to an increase 
in the size of the scattering units which in turn should give rise to a more 
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pronounced elliptic polarisation. Comparing the values of p, and Ap, 
given in Table I for the various glasses, one finds that the range of @ over 
which elliptic polarisation is detected, is directly correlated with Ap, and 
not with p,. It follows therefore that no reliable estimate of the size of the 
particles in two different scattering media could be had from a knowledge 
of the values of p, in the two cases. as they would be influenced to a large 
extent by the anisotropy of the particles. On the otherfhand, a comparison 
of the values of Ap, or the extent of ellipticity in the two cases would 
give us a better estimate of the relative sizes of the particles. Of the nine 
glasses investigated, No. | contains molecular aggregates of largest size and 
Nos. 4 and 8 contain those of smallest size. 


In a paper on the polarisation of light scattered by glasses Parthasarathy 
and others (1941) claim to have obtained negative results (i.e., p, was not 
less than 1) in a few specimens of glass some of which were highly strained. 
Their experimental method of determining the values of the depolarisation 
factors was not free from defects. It is well known that all glasses exhibit 
fluorescence. Fluorescence will enhance the value of p, as determined by 
photographic photometry. Their method for detecting the reversal of 
polarisation is far less sensitive than either the method of Krishnan (1938) 
or the new method of Sir C. V. Raman (1941) employed in the present 


investigation. Parthasarathy and others have worked with glasses having 
strains of macroscopic dimensions and tried to compare their results with 
Mueller’s theory [Hans Mueller (1938)] which on the other hand is based on 
the existence of frozen-in strains of microscopic dimensions. 


6. Conclusion 


The finite size in comparison with the wave-length of light, of the ele- 
ments of heterogeneity of any scattering medium is indicated by the existence 
of some anomalous depolarisation (i.e., 0, <1) and/or ellipticity in the scat- 
tering of polarised light in such medium. Of the two effects, the anomalous 
depolarisation or the reversal of polarisation is more easily observed and is 
very sensitive to the size of the scattering elements. Colloids such as gold 
sols, ferric hydroxide sol, vanadium pentoxide sol, etc., contain particles which 
are large enough to exhibit anomalous depolarisation but too small to show 
any measurable e'lipticity in the scattered light. The detection of anomalous 
depolarisation in any scattering medium is therefore a more delicate test 
for the existence of molecular aggregates than the detection of ellipticity. 
The one drawback for this test lies in the fact that the observed anomalous 
depolarisation may also be caused by the finite ang!e of convergence of the 
incident beam and/or secondary scattering. The effect of secondary 
scattering will be pronounced in the case of emulsion, protein solutions and 
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liquid mixtures very near the critical solution temperature. It is therefore 
necessary to eliminate the effects due to these disturbing factors in the deter. 
mination of p,. On the other hand, convergence of the incident beam and 
secondary scattering cannot give rise to any ellipticity in the scattering, as 
long as the medium is isotropic and contains particles which are small 
compared with the wave-length of light. Perrin (1942) has therefore rightly 
remarked that. “‘ the existence of some ellipticity in the scattering for an 
incident beam linearly obliquely polarised would be a much more sure test 
of the multipolar character of the scattering, and consequently of the 
non-negligible magnitude, in comparison with wave-length, of the elements 
of heterogeneity of the scattering medium ”’. 

The authors’ grateful thanks are due to Professor Sir C. V.. Raman for 
his kind interest in the work. 


Summary 


Sir C. V. Raman’s new technique for the study of light-scattering based 
on the use of a Babinet compensator, has been successfully applied for the 
detection of the anomalous depolarisation in optical glasses. Thus the 
existence of the reversal of polarisation when the incident beam is polarised 
with vibrations horizontal, has been definitely established by the new method 
which is of an objective character. The present investigation further 
reveals that in glasses the scattered light is elliptically polarised when the 
incident beam is linearly polarised in a direction oblique to the scattering 
plane. These observations support the view that molecular aggregates do 
exist in glasses. One important point which emerges out of the present 
investigation is that those glasses which give larger values of Ap,, exhibit 
elliptic polarisation to a correspondingly larger extent. Ap, is the difference 
between the observed value of p, and the anisotropic part of p, which can 


be considered to be equal to a = ; y It has been pointed out that for 


making comparative estimates of the sizes of particles in any two scattering 
media a knowledge of the values of Ap, in the two cases is by far more 
important than the corresponding values of p;. 
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Fics, | 6. Photographs showing anomalous depolarisation in glasses Nos. 6, 7, 4,9, 5 
and § respectively when the incident beam is polarised with vibrations horizontal, i.e., @ = 90°. 

Fics. 7 — 9. Photographs showing elliptic polarisation in the scattered light in glass No. 8 
for @ = 87°, 84° and 81” respectively. 
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1. Introduction 


It is well known that the principal vibration of the diamond lattice in which 
the two interpenetrating Bravais lattices of carbon atoms oscillate against 
each other manifests itself in the Raman spectrum as a sharp and intense 
line with a frequency shift of 1332 cm.-! at room temperature. Observations 
of the state of polarisation of this line were made by Ramaswamy (1930) 
who reported it to be strongly polarised. Bhagavantam (1930) in the course of 
his study of the Raman effect in several diamonds however, observed the 
line to exhibit only imperfect polarisation. This difference between the 
two observers is not surprising in view of the numerous disturbing factors 
capable of influencing the experimentally recorded result. More recently, 
Nayar (1942) made a re-determination of the depolarisation ratio, using a 
small octahedral crystal and immersing it in benzene to reduce the disturbances 
arising from the refractions and internal reflections at the surfaces of the crystal. 
He found that when the incident light was unpolarised, the two components 
of the scattered beam were of equal intensity, giving a value of 1-0 for the 
depolarisation factor. This, as he pointed out, differs from the theoretical 
value of 2 predicted by Saxena (1940) on the basis of Placzek’s theory for 
triply degenerate vibrations in cubic crystals. In view of these conflicting 
results reported by different observers and their divergence from the theo- 
retical values, it was considered desirable to obtain more accurate polarisa- 
tion data, taking the precautions necessary to eliminate the factors giving 
rise to errors, viz., birefringence in the diamond, refractions at the surfaces 
of entry and emergence, internal reflections and finally the convergence 
error present in all depolarisation measurements. The first three of these 
were avoided by choosing for the investigation from Sir C. V. Raman’s 
personal collection a diamond (N.C. 4) which had a symmetric and approxi- 
mately octahedral form and smooth lustrous faces and was free from any 
internal flaws or inclusions. The diamond was also free from any noticeable 
birefringence. This was tested by passing a beam of polarised light 
through the crystal and determining the state of polarisation of the emergent 
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beam. The disturbing effects of surface reflections and refractions were 
minimised by immersing the diamond in a high refractive index liquid, while 
convergence errors were reduced by using a small-aperture short-focus lens 
for condensing the light into the crystal. A careful study of the polarisa- 
tion characters of the Raman line as well as of the fluorescense bands of 
diamond has been made in the present investigation, using incident un- 
polarised as well as polarised light for excitation. 


2. Experimental Technique 


The diamond was kept immersed in a-bromonaphthalene contained in 
a small Wood’s tube just large enough to hold the crystal. The tube was 
provided with optically plane windows and was blackened throughout 
except for stops for the entrance and emergence of light, which besides 
cutting off parasitic light, served to confine the incident and scattered beam 
to mutually perpendicular directions. Light from a mercury pointolite arc 
was condensed by a small-aperture lens at the centre of the crystal through 
one of its octahedral faces, the scattered light being observed through another 
octahedral face approximately normal to the former. The scattered light 
was focussed on the slit of the spectrograph by a lens of small aperture. 
The convergence error in this arrangement was only 1%, the maximum 
semi-convergence angle being about I!°. A large aperture Ahrens Nicol 
interposed in the path of the incident light served as the polariser. A 
sodium nitrite filter isolated the 4358 radiation for exciting the Raman line 
effectively cutting off light of shorter wave-lengths which excite fluorescence, 
while during polarisation measurements of the fluorescence bands the filter 
was dispensed with. A quartz double-image prism of the Wollaston type 
properly oriented and placed just before the slit of the spectrograph split the 
beam into vertical and horizontal components. The spectra were photographed 
on Ilford Selochrome plates, the instrument used being a Hilger two-prism 
spectrograph having a dispersion of 28 A.U./mm. in the 4358 A.U. region. 


3. Polarisation Measurements of the Raman Line 


As usual in polarisation studies, three mutually perpendicular directions 
OX, OY and OZ in the crystal were chosen as reference axes; the direction 
of the incident beam was designated as along OX, the direction of the 
scattered beam as along OY while OZ was the direction mutually perpendi- 
cular to both. Since diamond is isotropic, it was possible to obtain only 
three different spectrograms of scattered light; viz.. with 

(a) incident light unpolarised, 

(b) incident light polarised with electric vector along OY, 

(c) incident light polarised with electric vector along OZ. 
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Quantitative measurements of the depolarisation values for these three cases 
were made and are given in Table I along with the theoretical values 
predicted by Saxena for the triply degenerate vibrations in the cubic class. 


TABLE [ 


E,?* 





Depolarisation 





Electric vector 


Direction of the Incident | 
| Experimental Theoretical 





AlongOY ... a i oad 1°5 
AlongOZ_ .. “e ar aie 0:6 


Unpolarised Light .. er 5 1-0 | 








| 








The depolarisation factor o was determined as usual by microphoto- 
metering the plates and obtaining the densities of the lines for the horizontal 
and vertical components. A set of intensity marks were obtained on the 
same photographic plate by the method of varying slit-widths, using as the 
source of illumination a small tungsten filament lamp fed from a battery 
at constant voltage. From the density-log intensity curve thus obtained 
for the wave-length of the Raman line, the intensity ratios of the horizontal 
to the vertical component were calculated. The ratios obtained in this 
way involve errors due to-the spectrograph caused by oblique refraction 
and reflection at the prism surfaces which strengthen the horizontal 
component. This error, which is characteristic of the instrument used, 
was determined by allowing light from a tungsten lamp after it had passed 
through a double-image prism to fall on the slit of the spectrograph and 
photographing the spectra for a series of exposures. The ratio of the inten- 
sities of the vertical and horizontal components was determined microphoto- 
metrically as described above. This gave a value of 0-71 for the 
spectrographic correction. The values in Table [ are corrected for the 
instrument error. 


4. Theoretical Derivation of the Polarisation Characteristics 
of the Raman Line 


Tables of selection rules for the different types of oscillations in crystals 
belonging to the seven classes of symmetry have been given by Placzek (1934, 
1938). On the basis of these, Saxena (1940) has calculated the polari- 
sation ratios and the total intensities of Raman lines in terms of the tensor 
components in the transverse and longitudinal scattering. He has shown 
that Raman lines in crystals may be considered as due to the deformations 
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of the optical polarisability ellipsoid of the solid produced by the internal 
oscillations and that the deformations can be expressed in terms of the 
six components of the change of polarisability tensor. The equation for the 
deformed ellipsoid is given by 


(A+ €22) X2+ (B+ €,,) yP (CH dg) 29+ ey XY Ezy YZT Eqs XZ= I 
where A, B and C are the lengths of the principal axes of the polarisability 
ellipsoid and €;; = (€xx» €yy» €:2) €i&= (Exy» €ycr x) are the six Components 
of the change of polarisability tensor. It is possible to evaluate these tensor 
components by illuminating the crystal along a fixed direction, say OX, and 
observing the scattered light along OY, and calculating the induced moments 
in the scattered radiation. The components of the scattered radiation lie 


along the x and z directions and the induced moments in these directions are 
given by 


AM, = «,° E, + ©, E, +«° E, 
AM, =«,, - E, +«4,° E, + «,° E, 
Since the depolarisation ratio p is given by 
_ (AM,)? 
Pp (AM. 
it will be seen that for an incident electric vector along OY 


_(éxy)* 

Poy “a <,,)8 

and for an incident electric vector along OZ 

2 

roe (i 

while for incident unpolarised light 

Puno: = (<,»)* + (<,.)*. 
_ (<,,)° + (€..)* 

The principal vibration of the diamond lattice falls under the triply 

degenerate class for which the selection rules given by Placzek are 
Exe = Eyy = Exe = 0 
€xs = €ys = Exy 0 

Thus we find that poy= 1, poz = oo and p,,,.= 2. 

It will be seen from Table I that while the results for incident 
unpolarised light agree with those obtained by Nayar (in spite of the small 
error introduced in his results by the crystal birefringence), they agree 
in no case with the theoretical values obtained on the basis of the selection 
rules given by Placzek. This discrepancy between the theoretical and the 
experimental values indicates the necessity for a reconsideration of the theo- 
retical ideas, on the basis of which the selection rules have been derived. 
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5. Polarisation of the Fluorescence Bands of Diamond: 


Table Il gives the values for the depolarisation ratio p for the 
fluorescence bands of the 4152 system for the 3 cases mentioned above, 
corrected for the instrumental error. 


‘ 


TABLE Il 








Incident light unpolarised .. 8 s = am. te 1-0 








Incident light polarised with electric vector along OY .. i ne 0-8 





Incident light polarised with electric vector along OZ .. ae ne 0-2 








Ramaswamy had previously observed that with incident unpolarised light 
the fluorescence bands were depolarised. It will be observed from the table 
that the two components of the fluorescent beam are of equal intensity for 
incident unpolarised light, while the bands appear partially polarised with 
the incident electric vector along OY and strongly polarised with the incident 
electric vector along OZ. No obvious difference was noticed in the polari- 
sation characters of the main band at 4156 A.U. and the subsidiary bands at 
longer wave-lengths. 


In conclusion, the auther wishes to express her grateful thanks to Professor 
Sir C. V. Raman for his helpful interest and encouragement in this work. 


6. Summary 


The depolarisation values for the principal Raman line in diamond 
have been determined using a natural crystal and taking the necessary 
precautions to eliminate errors due to the birefringence of diamond and to 
refractions and internal reflections in the crystal. The results obtained for 
the three different cases, viz., with the incident light unpolarised and polaris- 
ed with the electric vector along the OY and OZ axes. are found to differ 
notably from the theoretical values predicted by Saxena on the basis of the 
selection rules given by Placzek. 


The polarisation of the fluorescence bands of the 4152 system has 
also been studied. It is found that with incident unpolarised light the 
bands appear depolarised while they appear polarised in the other two cases. 
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NORMAL OSCILLATIONS OF THE T, CLASS 
- DIAMOND STRUCTURE 


By S. BHAGAVANTAM 


(From the Department of Physics, Andhra University, Guntur) 
Received July 25, 1944 
1. Introduction 


IN an earlier paper,' group theoretical methods have been used by the author 
for obtaining the normal oscillations of the O, class diamond structure and 
the results appropriate to a repeating unit which has eight times the volume 
of the smallest unit cell were given. Such a repeating unit is formed by 
taking twice the primitive translation as the side of the cell along each direc- 
tion instead of the primitive translation itself and contains 16 carbon atoms 
whereas the smallest one contains only 2. 


It has been recently suggested by Raman? that the commoner type of 
diamond has only tetrahedral symmetry and should be placed under the T, 
class. It is considered desirable to work out the normal oscillations 
of this structure also. In the present paper, results of applying group theo- 
retical methods to such a structure are given. So that the results may be 
applicable to other tetrahedra! crystals also (T, class), a diamond arrange- 
ment with half its atoms different from the rest is assumed. Such an assump- 
tion automatically secures the absence of a centre of symmetry. The 
repeating unit is again formed by taking twice the primitive translation as 
the side of the cell along each direction instead of the primitive translation 
itself. Such a cell contains 8 atoms of type I and 8 atoms of type II 
whereas the smallest one contains only one of each. 


2. T, Class Diamond Structure and its Symmetry Operations 


In Fig. | is shown a portion of this structure which is made up of two 
interpenetrating Jattices. The dark circles denote type I atoms whereas 
the white circles denote the type I] atoms. The smallest unit cell is a rhom- 
bohedron formed by the primitive translations 1,2; 1,3 and 1,4. Such a 
rhombohedron contains only two distinct atoms which are numbered | and 5. 
In such a case translations 1,2; 1,3 and 1,4 cannot be regarded as distinct 
from the identity operation and all the atoms numbered 2, 3, 4, 9, 10, 11 and 
122 
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12 are equivalent to 1, whereas those numbered 6, 7, 8, 13, 14, 15 and 16 are 
equivalent to 5. The symmetry group will then consist of 24 elements only 
and may easily be dealt with.* The unit cell chosen in this paper is, however, 
a larger rhombohedron formed by the primitive translations which are twice 
1,2; 1,3 and 1,4. Such a rhombohedron contains 16 distinct atoms which are 
numbered | to 16 in Fig. 1. Besides the 24 elements of the simple group, 
seven translations (1,2; 1,3; 1,4; 1,9; 1,10; 1,11; 1,12). which were hitherto 
identical with the identity element, have now to be regarded as distinct 
symmetry operations. These, along with the identity operation constitute a 
sub-group of order 8 and the total group of symmetry operations appropriate 
to this structure in which there are 16 non-equivalent points as shown in 
. Fig. 1, will be formed by obtaining the product of the simple 24 elements 
with the above translational group consisting of 8 elements. The resulting 
group is of order 192. The elements of this group fall into 13 conjugate 
classest and the appropriate character table is given here. The notation 
used is similar to that employed in the earlier papers of the author. The 
table shows that besides the translation, there are three threefold, two six- 
fold, two fourfold and two eightfold degenerate normal oscillations coming 
under various symmetry classes. Translation comes under the class F,. 


If all the last eight columns and the last eight rows are deleted, we get 
the character table containing only five conjugate classes and if the number 
of elements in each conjugate class is properly adjusted, we get results that 
are appropriate to the simple group of 24 elements which refers to the case 
where the smallest cubic cell containing 2 atoms only is regarded as the 
repeating unit. 


> 


3. Normal Modes and Normai Frequencies 


Below is given one representative mode under each class. Those omit- 
ted may be written down from considerations of symmetry and taking into 
account the order of degeneracy in each case and the manner in which each 
normal co-ordinate transforms under the different symmetry operations. 
If the normal co-ordinates given below are read with reference to Fig. 1, it 
is easy to get a physical picture of each one of the modes. 


* Identity, eight trigonal axes, three digonal axes, six planes of reflection and six tetragonal 
rotation reflection axes constitute these 24 elements. 
+ Writing down all the 192 elements in the form of circular permutations and classifying 


them into conjugate classes is a very laborious process but follows the well-known methods of 
group theory. Details are not given here as they would occupy much space. 
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Mg (X;, + Xt Xgt Xq+ Xt Xyo+ Xia + X12) 

— My (Xs Xgt X_pt Xyt Xigt Xyat X53 + X16) pe F, 

(x, + X—)— (X7 + Xe) + ist X14)— (X15 + X26) s - F, 

(x, + X2)— (X%3+ X4) + (X%q+ X49) — (X11 + X12) “+ . F; 
(X4— Xgt+ Xy— Xq) t+ (Xy—Xiy t+ X11— Xia) -- 
(¥s— Ye + Ya— Ya) + Nis— Via t Vis— Vis) -- 
(p2+ Pst Pat Ps)— (Pit Prot Put Pr) -- 
(Po+ Pz + Pat Pis)— (Ps+ Pris t Pis+ Pre): - 

~ (X,— 21) + (X2— Ze+ Xg3— Zg+ Xq— 24) + (X_— 25) 

— (Xy9— 219+ Xa — 211 + X12—212) 


H, 


K, 


(x;— 25)— (Xg— Zp + X_— 27+ Xg— Zg)— (X13— 213) 
+ (%1a— 214+ X15— 215 + Xis— 216) °° ++] 


| ow 


p in K, stands for a displacement x, y, z. my, and mz, are respectively the 
masses of the type I and the type Il atoms. The frequencies of these modes 
can be evaluated in the usual manner. Below is given a statement which 
shows the frequency in each case and contains a description of the corres- 
ponding normal mode.* A stands for 47’v?. yw is the reduced mass and 
i = «at uf K, is the force constant referring to any pair of atoms which 
constitute the nearest neighbours. K,’ and K,” are the force constants 
referring respectively to the nearest pair of type I atoms and of type II atoms. 
K,’ and K,” respectively represent the forces called into play when the angle 
between any two valence bonds which meet at the type I and the type II atom 
varies. p stands for the distance between a pair of nearest neighbours and 
is the same as the length of the valence bond. 


It may be noted that if m, and m, are both replaced by m, the mass of a 
carbon atom, K,’ and K,” are both replaced by K, and K,’ and K,” are 
both replaced by K;, the structure reduces to that of O, diamond as a centre 
of symmetry is automatically introduced. In sucha case it is easily veri- 
fied that the fequency of F, becomes identical with that of F. in O,; diamond.! 
The frequencies F; and F; become identical and. they together correspond 
to the sixfold frequency H, of O, diamond. The pair coming under H, 
pass over into the two coming separately under H, and H, in O, diamond. 
Similarly, the pairs coming under K, and M, pass over respectively into those 





* K and p occurring here have a different significance and are not to be confused with those 
occurring in the normal modes. 
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coming separately under K, and K, and M, and M, in O, diamond. One 
distinguishing feature of the T, structure is the activity of the principal 
oscillation coming under the F. class both in Raman effect ahd infra-red 
absorption. 

4. Summary 


Group theoretical methods have been applied for obtaining the normal 
oscillations of the T, class diamond structure on the basis of a 16 atom cell 
as the repeating unit. It is shown that besides the translation, there are nine 
normal oscillations three of which are threefold degenerate, two of which are 
fourfold degenerate. two of which are sixfold degenerate and two of which 
are eightfold degenerate. Expressions are derived for the corresponding 
normal frequencies by postulating suitable potential energy functions. 
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THE HYDROGEN BOND AND DIAMAGNETISM 


By S. V. ANANTAKRISHNAN AND P. S. VARADACHARI 
Received April 21, 1944 


1. Introduction 


THE interpretation of susceptibility data in terms of atomic contributions 
with a few “constitutive corrections” developed by Pascal'® has been so 
successful that even recent studies make no attempt to depart from this rule. 
Pauling’ is constrained to remark ‘* Diamagnetic polarisation can be cal- 
culated with reasonable accuracy as the sum of terms corresponding to the 
atoms composing the substance with, perhaps, some small corrections for 
the type of bonds involved. Diamagnetic polarisation has not been found 
to be of much value in the investigation of the structure of molecules.” How- 
ever, attempts have been made by Gray and co-workers! ™ to interpret dia- 
magnetic data in terms structure and though * calculated ’ and * observed’ 
values show remarkable agreement, the coincidence as will be shown in a 
jater communication, is largely fortuitous. 


Among structures that influence physical properties of organic com- 
pounds, hydrogen bonds have been known to cause noticeable deviations 
from usual laws while the methylene group has been found to exert a remark- 
ably constant influence on these properties when one compares the mem- 
bers of different homologous series.??° One may reasonably expect, 
therefore, that, taking these two factors, the extent to which hydrogen bond- 
ing alters diamagnetic susceptibility may be gauged. Other investigations” 
show that hydrogen bonding may be intramolecular or intermolecular while 
a weaker type has been designated by Bernal and Fowler® the * hydroxyl 
bond ’. 


Apart from a preliminary report by Angus and Hill! there does not 
appear to be any special study of the problem and the present survey is 4 
prelude to a systematic structural interpretation of diamagnetism. In the 
comparison of susceptibility values, the question of choice between the 
results of different workers arises. For a proper comparative study, @ 
necessary and sufficient condition is that values should be self-consistent 
and any deviation to be interpreted must be bevond the limits of experi 
menta! error. These limitations preclude any correct interpretation in 
some cases even if other physical properties indicate the existence of 4 
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hydrogen bond. The results of Angus and Hill (/oc. cit.), the only magne- 
tic investigation of the hydrogen bond, illustrate this point as may be notic- 
ed from their data presented in Table I for benzoic and salicylic acids. 


TABLE I 








Benzoic Acid Benzene : 0-702, 0-720 | 0-702 0-004 
. 0-7020 0-685? 


0-698, | 0-650 


| 
Substance | Solvent c (wt. %)| Solution Solute Solvent 
{ 








Ethyl acetate . 0-606, 0-432 0-612; 0-003 
: 0-6050 0-483 


| 
Salicylic Acid : ; Chloroform ° 0-516; 0-632 0-515 0-002 
° 0°5160 | 0-582 


| Ethyl acetate | 4. | 0-6120 0-602 
| ial 














0-6077 | 0543 





| 





[t will be noticed that all the susceptibility values are for very dilute 
solutions and lie within the limits of possible experimental error. From a 
later section it will be found that the change in diamagnetism of a polyatomic 
molecule through hydrogen bonding cannot be as large as one might infer 
from the x,;u. Values and the variations with even changes of concentra- 
tion appear to be abnormally large. [It may be mentioned, however, that the 
changes observed in the case of solutions of benzoic acid are in the right 
direction. Judging from these results. one has to conclude that the study 
of dilute solutions may not be a fruitful line of approach to the problem: of 
hydrogen bonding. A different approach is therefore indicated. 


2. Intramolecular Hydrogen Bonds 


Well-known examples of intramolecular hydrogen bonding are found 
in aromatic compounds. Taking closely related structures differing only by 
the presence of a hydroxyl group, molar susceptibility values may be expected 
to indicate hydrogen bonding. Though phenol is considered to be associated, 
the work of Badger and co-workers**® indicates that no serious error will be 
introduced by taking it as unassociated for our present purpose. The sus- 
ceptibility data for a few closely related structures are given in Table II. 


The first two pairs of compounds indicate that the replacement of a 
hydrogen by hydroxyl in a six-membered ring brings about the same change in 
susceptibility. Substances A in the latter pair are known to involve hydrogen 
bonded structures and the extra contribution to diamagnetic susceptibility 
may naturally be traced to this factor. The net increase in the value works 
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TABLE II 
(Values taken from I.C.T.) 











Substance A Substance B X, Xx. Xa, | ae %u, Xn, 
— Sais : 
Cyclohexanol Cyclohexane | Q+-726 0-810 72+6 | 67-3 53 
Phenol Benzene 0-648 0-712 60-9 | 55-4 5-5 
Methy! salicylate Methy! benzoate 0-580 0-602 88-2 | 81-9 63 
2: 4-Dinitro phenol | m-Dinitro benzene} 0-397 0-398 73-0 | 67-0 6-0 




















out to nearly one unit per mole and incidentally indicates this as an alternative 
method that can be advantageously used in the study of hydrogen bonding 
especially in compounds of low molecular weights. 


3. Intermolecular Hydrogen Bonds 


If, as has been noticed in the previous section, hydrogen bonding in- 
creases the diamagnetism, the reverse process of rupture of hydrogen bonds 
should lead to a decrease. The conversion of an aliphatic carboxylic acid 
to the corresponding methyl ester involves not only the addition of a methylene 
group but also the rupture of the hydrogen bonds in the dimer found normally 
in all liquid and solid aliphatic acids at room temperature.'® Consequently, 
the apparent contribution to the molar susceptibility by the methylene group 
should appear to be smaller than the normal value by about one unit as 
the change involves the rupture of one hydrogen bond per molecule. 


Ignoring the very low value reported by Bhatnagar and co-workers’ 
it is generally found that the increase in y,, caused by the addition of a — CH, 
group is 11-69 x 10-* per mole and one may justifiably assume the constancy 
of this value for calculations. The agreement among the different workers 
may be seen by examining the data in Table III and the curves in Fig. 1. 




















TABLE III 

| az: 

XAng x X1.6,7, Xan, | x x 

Acid Qy | «ay | Ge | (33) used | 
| L cacaidaalal 

Formic ” = me fee 0-432 0-433 0-433 19-88 
Acetic - || 0+528, -. | 0526 | «(0-525 0-528 | 31-72 
Propionic S | 0+585, | | +587 ~ 0-585 | 43-36 
Butyric eo ..| 0-626, | 0-625, | 0-632 .. | 0627 | 55-20 
Valeric ie Baer 0-654, |... oe 0-655 | 66-85 
Caproic ¥ ..| 0°676, | 0-673, | .. mi 0-676 | 78+14 
Heptylic s — es 0-600 | «| . 0-690 | 89-74 
Caprylic m spo | 0-7053 } | . 0-705 | 101-6 








Plotting X, against the number of methylene groups n (CH,) it will be 
noticed from Fig. | that all the points lie on a straight line and by extrapola- 
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100 






































’ ' “No of CH, groups , 

Fic. | 
tion from the curve, the diamagnetic contribution per methylene group works 
out to 11°69. Substitution of hydrogen by CH, in the alkyl part of the 
carboxylic acid does not, therefore, imply any abnormality. It will be noticed 
from Table [V and Fig. 1 (a) that the replacement of carboxyl hydrogen is 
different. 











TABLE IV 
aes | | Difference 
| aoe 
ja x,, eSter | x, ester | calculated 
Acid : “ie calculated | iia siead observedt and observed 
_ | | %y values 
| | | of esters 
—_————$— Sn a —— _ — 
Formic | 19-86 | 31-52 (0+5267) 31-60 | —0-08 
Acetic we el 31-72 | 43-36 0-5719 42-34 | 0-98 
Propionic ste en 43-35 | 55+04 0-6136 | 54-04 1-01 
Butyric at oil 55+20 66°84 | 0-6446 | 65:79 | 1-05 
| 


| | { 





* x,, is calculated assuming the monomeric form of the acid alone. 
+ The value for formic ester is that given by Rao and Narayanaswami.** Higher values 
reported can be traced to water present as impurity which is difficult to remove in this case. 
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Angus (Jvc. cit.) has observed that in changing from the acid to the 
methyl ester, the diamagnetic contribution of the methylene group is 10-69 
or one unit less than the normal value. 


Assuming the constancy of the methylene group contribution, it will 
be noticed from the last column that the difference between the calculated 
and observed values of the esters is just one unit per mole, a difference that 
may be expected for a dimer with one hydrogen bond per molecule (vide 
Section 6). The first member behaves abnormally. Like all other carboxylic 
acids, the dimeric form of formic acid has structure I, the C-H hydrogen 
being free.* If the acid is converted into the methyl ester, the hydrogen bonds 
of the carboxyl group is destroyed but.a new set of dimers having the alterna- 
tive structures IT or IIT involving the C-H hydrogen is presumably formed.+ 


CH; 
l 
> H-O. >H H 
p> \ yo o> 
H-C C¢-H O=C C=O H,CcO-C C-OCHs 
No-H «0% \H 0% \u <—07 
| 
CH, 


(I) (II) (111) 


Since the number and nature of the hydrogen bonds apparently remain 
the same in both the acid and the ester, the contribution by the methylene 
group appears normal. Methyl formate is thus abnormal when compared 
with the other methyl esters and the point corresponding to methyl formate 
on the X,—nCH, curve (Fig. 1a) does not fall in line with the rest as it 
happens to be dimeric involving a hydrogen bond. | 


In a later section it will be shown that the provable structure is II and 
not III since it involves C=: O as in the acids and as such may be expected 
to have the same effect on susceptibility. The Raman spectrum with a modi- 
fied C-H frequency'® *4 confirms this view and the carbony! frequency also 
supports the same idea. 


Mention should also be made here of Angus’s observation? that the 
change in susceptibility noticed between formic acid and the formate ion 
is due to hydrogen bonding in the acid. From our foregoing analysis it 
will be noticed that the change caused by hydrogen bonding in the molar 


* The Raman spectrum of formic acid*® shows a normal C-H frequency which will not be 
the case if it were involved in hydrogen bonding as in the complex of the type postulated by 
Murti and Seshadri.'* The feeble second carbonyl frequency is traceable to the known exist- 
ence of about 10% monomer even at room temperature.}* 


+t Methyl formate has higher m.p. than methyl acetate indicating complexity. 
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susceptibility is small though definite. The large observed difference be- 
tween the acid and the ion has probably arisen more from the structural 
changes in the carboxyl group. The non-equivalence of the resonance struc- 
tures and presumably different carbon oxygen bond-lengths in the case of 
formic acid and the equivalence in the formate ion may give rise to the 
observed difference. This view is supported by the Raman spectrum which 
indicates no carbony! frequency in the case of formates.*® 


4. Alcohols and Water 


Unlike the cases considered in the previous sections, the alcohols and 
water (which may be considered the first member in the homologous series 
R-OH) show variations of a different type. It is recognised that all abnor- 
mal physical properties of this class of compounds result from association 
through hydrogen bonding.?? In order to differentiate these weak polarisa- 
tion bonds, Bernal and Fowler® have in their theory of a quasi-crystalline 
structure of liquid water postulated a hydroxyl bond. The work of 
Cabrera and Fahlenbrach® on water and a series of alcohols as well as that 
of Auer on water® seen to indicate a small but definite increase in diamag- 
netic susceptibility with temperature (cf. also Varadachari®*). Since rise 
of temperature can reasonably be expected to favour depolymerisation, 
this implies that the monomer has a higher diamagnetic susceptibility than 
the dimer, a result apparently the opposite of what has been noticed in the 
case of carboxyl compounds. This will be discussed in a later section. 


5. Inorganic Oxy Acids 


If one assumes as in the organic acids that the methylene group con- 
tribution to diamagnetic susceptibility is a constant. the observed suscepti- 
bility of dimethy! sulphate (x,, 62-12)!* leads to a X,, value for sulphuric acid 
of 38-84. The careful investigations of Farquharson" gives an experimental 
value of 42-3 for pure sulphuric acid confirming the value given in Inter- 
national Critical Tables. This large difference may be attributed to hydrogen 
bonding and other physical properties also indicate the possibility of larger 
aggregates in 100% sulphuric acid. 

The case of iodic acid is even more interesting. The x, value changes 
from 47-01 in the crystalline orthorhombic form to 41-47 in solution. This 
large difference in the susceptibility is partly due to the break up of the 
hydrogen bonds present in the crystal and partly to ionisation in solution. 
The investigations of Helmholtz and Rogers'® have shown that the crystal 


with a unit cell of four iodic acid molecules contains bifurcated hydrogen 
bonds 
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A proper discussion of the contribution of hydrogen bonds to the dia- 
magnetism of inorganic acids is, however, not feasible at present especially 
since most of the data are derived from the measurements with correspond- 
ing salts and as such unsuitable for interpretations of this nature. 


6. The Nature of the Hydrogen Bond and its Contribution 
to Diamaznetism 


Since diamagnetism is liable to be influenced by the interaction of 
orbital electrons, the change in diamagnetism in hydrogen bonded structures 
should be related to the nature of such bonding. Earlier views implied 
an apparent increase in the covalency of hydrogen but this has been severely 
criticised. The hydrogen bond may now be treated as an essentially hybrid 
structure, the resonance involving hydrogen attached to one atom by a true 
covalency and the other by an electrovalency. 


In any structure 
A—H....B 


the distance between the atoms A, B taking part in the hydrogen bond is 
invariably found to be less than in the case for an isolated pair of atoms and 
the shorter this distance the greater the strength of the hydrogen bond. The 
structure may be more correctly pictured as essentially arising from the loca- 
tion of a proton at an equilibrium distance with the bonding electrons more 
spread out between the two extreme positions. 


Following Van Vleck?’ the diamagnetism of a polyatomic molecule may 
considered to involve a paramagnetic term arising from the distortion of the 
orbital electrons and a diamagnetic term. Both these are independent 
of temperature. The formation of a hydrogen bond can influence the dia- 
magnetic susceptibility in four directions :— 


1. Increased diamagnetism through the sperading of the electron cloud 
or the formation of larger molecular orbitals of the bonding electron. 


2. Decrease in the paramagnetism associated with the bonds through 
a release of strain in the groups involved. This is likely to be appreciable 
only in the cases where multiple linkages are involved. 

3. The formation of an additional bond make a small paramagnetic 
contribution through the distortion of the electron orbits. 

4. A large diminution in diamagnetism resulting from a distortion of 
the orbits and a decrease in the effective area through the presence of 4 
proton inside the electron cloud of two atoms connected by hydrog2n 
bonding. 
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Any change noticed that has to be a resultant of all these factors. 


An analysis of the structures involved in hydrogen bonding shows 
that we can broadly divide them into two categories according as multiple 
bonds or single bonds are involved. Examples of the first type are provided 
by carboxylic acids, esters of salicylic acid, nitrophenols and sulphuric acid, 
all of which involve structures of the type 

O—H....0=A 
where A stands for carbon, nitrogen or sulphur. As mentioned earlier, 
the dimeric forms of the carboxylic acids are usually assigned the structure 


f°" iat a 
Cc 2 


No- ane of 
As a consequence of hydrogen bonding, the bond strength of both C =O 
and O-H are affected and this is noticed clearly from the shifts in the car- 
bonyl frequencies of Raman spectra!* 2° and the carbonyl and hydroxyl 
frequencies in the infra-red studies of Badger and co-workers.® 


From the large paramagnetic contribution attached to a double bond, 
any diminution in the double-bonded character (a decrease in strength of 
binding) due to hydrogen bonding must lead to a decrease in the Van Vleck 
paramagnetic term. The observed net increase in diamagnetism is then 
a clear indication that the dominant factor in this case is presumably the 
diminished paramagnetism. One can also visualise the ring structure as 
involving a mobile electron system with a continuously changing bond value 
throughout the ring. Just as in the case of aromatic ring systems this may 
result in an increased diamagnetism.'® The study of anisotropy of the 
crystalline!” acids may be expected to throw further light on this aspect of 
problem. The intramolecular ring structures considered in the earlier sec- 
tion also belong essentially to the same type, a consequence of which is seen 
in the same molar increase in diamagnetism of one unit per mole. per 
hydrogen bond. 


The principal change in the case of iodic acid from the solid to solution 
is traceable to a break up of the bifurcated hydrogen bond. This special 
structure which may be taken to imply a larger electron orbital would 
naturally lead to an abnormal increase in diamagnetism. It is also possible 
that the multiplanar system in the unit cell of four molecules constitutes a 
mobile electron system with the consequent implications. 

Contrasted with all these structures, we have the alcohols and water 
in which the break up of the hydrogen bond leads to an increase in dia- 
magnetism. Structurally the polymers involve only single bonds: 
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It is generally recognised that every new bond increases the distortion of the 
electron system, thereby increasing the paramagnetic term while the dimi- 
nution in the paramagmetic term (factor 2) present in the acids is absent. 
Possibly, in the case of water and the alcohols the increased paramagnetic 
term outweighs the diamagnetic term associated with the increased orbitals 
of the hydrogen bonded system. The apparent difference in behaviour 
between the two types then becomes understandable. 


Interpreted thus, we have to admit that any observed change in suscep- 
tibility is determined primarily by the groups involved. This may then be 
of value in deciding between alternative structures as in the case of methyl 
formate (vide Section 3). From the observed diamagnetic susceptibility, 
structure II given in that section appears most likely. 


7. Summary 


A critical analysis of available data shows the constancy of magnetic 
contribution by the methylene group. Assuming this value it is shown that 
hydrogen bonding leads to increased diamagnetism, the susceptibility change 
being one unit per mol whenever structures of the type O—H....0=A 
are involved. A bifurcated hydrogen bond as in iodic acid apparently leads to 
an even larger value but a proper evaluation of this is at present complicat- 
ed by other factors. The case of alcohols and water which show a decrease 
in diamagnetism on association indicates that there is here an increase in 
the paramagnetic term. associated with distortion. The use of this in inter- 
preting structures is indicated. 
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1. Introduction 


IN a series of papers published in these Proceedings (Dayal, 1944 a-f) the 
author has successfully evaluated the specific heats of simple solids, and 
the thermal expansion of three common alkali halides, :iz., NaCl, KCI and 
KBr, on the basis of the Raman dynamics of crystal vibrations. Of the 
latter class of compounds, the fluorides of lithium and sodium are particularly 
interesting. They are constituted of lighter ions and as such, are expected 
to have higher characteristic frequencies of vibration. The deviations from 
the Debye specific heat formula, noticed in the specific heat curves of the 
other alkali halides and of most simple solids will, therefore, in their cases 
occur at higher temperatures. Measurements of the specific heat and 
of the thermal expansion are only available for lithium fluoride. Its specific 
heats have been measured by Goldmann at Gottingen but the results do 
not seem to have been published. Blackman (1935), who had access to 
these results, has pointed out that, with the lowering of the temperature, 
the Debye characteristic temperature at first drops to 610° K., but later rises 
sharply, at still lower temperatures, to 690° K. without showing any tendency 
to reach a constant value. The coefficients of thermal expansion have been 
measured by Adenstedt (1936) at low temperatures, and by Eucken and 
Dannohi (1934) at higher ones. The former, who has quoted a few of 
Goldmann’s measurements of specific heats, has found that the ratio of the 
thermal! expansion to the specific heat decreases continuously at the low 
temperatures. 


In the present paper, it is proposed to evaluate the specific heats and 
the coefficients of thermal expansion of these fluorides on the basis of the 
Raman theory of crystal vibrations. In the case of lithium fluoride, a satis- 
factory agreement is shown to exist between the theoretical and the experi- 
mental values of the specific heat and the thermal expansion. 
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2. Evaluation of the Frequencies 


In the paper on the specific heats of the alkali halides (Dayal, d), the 
repulsive force constant ¢, was calculated from the compressibility or the 
infrared absorption frequency of Barnes, on the assumption that the non- 
coulomb forces between the atoms other than the neighbours are small and 
can be neglected. In the case of rock-salt, its value was taken from 
Lyddane and Herzfeld’s paper (1938), where it was evaluated from Born 
and Mayer’s potential (1932). Practically the same value of c, can, however, 
be obtained, either from the infra-red frequency, or from Durand’s value 
of the compressibility (3-8 x 10-'* cm.?/dyne) at low temperatures. In 
the case of lithium fluoride, the compressibility is uncertain, having been 
given as 1-17 x 10-12 by Bridgman (1932) and 1-53 x 10-1® by Slater (1924). 
If we use the infra-red frequency to obtain this constant, there is a small 
but systematic deviation between the calculated and the observed values 
of the specific heats. Obviously, the above approximation is not enough 
and it is necessary to include the effect of the forces between the next nearest 
neighbours as well. We also expect the presence of non-central forces from 
its resemblance to the alkaline earth fluorides in its sparing solubility, and 
from the fact that the sum of the two ionic radii, as given by Goldschmidt 
(1929), is slightly greater than the lattice constant. The presence of weak 
forces of this type must also be presumed in the case of the other alkali 
halides because their elastic constants fail to obey the Cauchy relation at 
the low temperatures. If the elastic constants of both the fluorides were 
known, the contribution of these to the force constants could be directly 
calculated from them. In the absence of this information, we shall ignore 
their existence and confine ourselves only to the central forces between the 
neighbours and the next nearest neighbours. Following Born and Mayer 
(l.c.) and later workers, we assume that the repulsive potential between the 
atoms is given-by an exponential function of the type be~’’’, where + is the 
distance between the atoms and b and p are constants. 4 is usually obtained 
from the ionic radii, but the constant so obtained does not always give the 
correct value of the compressibility. For example, Huggins (1937) calcu- 
lated the compressibility of lithium fluoride in this manner and obtained a 
value -88 x 10-12, which is quite different from either of the experimenta: 
values quoted above. We will, therefore, make no assumption as regards 
the values of the constants b and ¢. Since the repulsions between the metallic 
ions are very small on any theory, we neglect them and consider the repul- 
sions between the fluorine ions only. The Van-der-Waal’s forces are also 
neglected. The potential energy ¢ per cell can be written as 

2 


BK nie =e + 6v (ro) + 6v* (ro), (1) 
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where a is the Madelung constant, e the ionic charge and ry, the distance 
between the nearest ions at the absolute zero. u(r) is the repulsive potential 
between the nearest alkali and the fluorine ions and is of the form 
b exp (— rip). v* (r) denotes the repulsion between the nearest fluorine ions 
and is of the im b exp(— rg/2/9’). The following definitions are 
introduced for convenience. 

= — d2v “) 

ro wr" Jr=r 


sa 2 ro *(* -) 


4rre* ail d2vF (r) Q) 


To? 








4ne® 2 sdv¥ (r) 
ye” ( wr Jr=rnVv2 

The four force constants ¢,, etc., are subject to the following three condi. 
tions. In the first place, the equilibrium condition of the lattice must be 
df 
or 
from the observed values of the infra-red frequency and the compressibility. 
These are given respectively by the formule 


satisfied, i.e., must vanish. The other two conditions are obtained 


> >  4ne? ] 
4r* yp" a= ure? [4 TT 2¢5 —~ 6 | (3) 
l N 20d 4 dd 
and 3 OV [ a | (4) 


where e is the ionic charge, # is the reduced mass, 3, the compressibility and 
N and V stand for the Avagadro number and the molecular volume res- 
pectively. 

The four force constants are reduced to three if we neglect c., which is 
a very small quantity, not more than one per cent. of c,. The three condi- 
tions are now enough to determine them completely. If we take the average 
of Bridgman and Slater’s values of compressibilities for LiF and allow for 
about ten per cent. change from room temperature we. get its value as 
1-25 x 10-2 cm.?/dyne at the absolute zero. With this value of the 
compressibility and taking the infra-red frequency of Barnes (1932), all the 


conditions are satisfied if * is nearly equal to 7:6. We have actually 


assumed this value for ry 9 and given a small value to c,, which is just enough 
to satisfy the equilibrium conditions. For NaF, the observed compressibility 
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is 2:11 x 10-1? cm.*/dyne (Bridgman, /.c.). Allowing for about fifteen 
per cent. change, 1-75 x 10°!” is a reasonable value at the absolute zero. 
With this value of the compressibility and taking v¥ =0, and r/p =9, 
all the above conditions are fulfilled for this salt. The method of calcu- 
lation of the frequencies, therefore, reduces to the same form as was dis- 
cussed in the earlier paper (d). For LiF, the numerical values of the con- 
stants are c, = °337, cy = — :0437, cy, = -026 and c,, = — -002. In 
the case of NaF, c, = -421, c. = — -0464. The room temperature values 
of ry» have been used throughout the calculations. The repulsive parts of the 
coupling coefficients as defined by Kellermann (1940), have been recalculated 
for all the frequencies of the Raman theory. The new values are 

R i * e 4re? 


SE : Cc Co + (Cyp + Cop) X 
xx ro? [c, + 2 (Cyr or) 


x (2 — cos ¢, cos d, — cos 4, cos 43) + 2Csy (1 — cos d, COS d5)] 
11 4re? ; , 
7 | a | = ctl [(Ciy — Cog) Sin , sin dg] (5) 
xy Ve 
where ¢, = 2k,rg 4, =2tkyry $3 = 27k,7o. 

The coefficients involving different types of atoms (1, 2) and those 
involving the metallic ions only, remain the same as before. The coefficients 
given in (5) are substantially’ the same as calculated by Thirring (1914), 
except for the small quantities c, and c.,. which, for obvious reasons, do not 


find a place in Thirring’s equations. k,, k,, k, are the components of the 
> 


vector a A being the wave-iength of vibration, and are identical with 


Tx» Ty» T;, used in the paper on the specific heat of the alkali halides (4). 
The frequencies have been calculated by the help of the expressions given 
by Kellermann and are given in Table I. The electrostatic coupling co- 
efficients have been taken from the same work. The frequency designation 
is the same as in the earlier papers. For the convenience of reference, the 
name of the vibrating ion in the vibrations of the octahedral plane has also 
been given. The other ion remains at rest in this mode. 


3. Evaluation of the Atomic Heats and the Coefficients of 
Thermal Expansion 


The molecular heats calculated from the frequencies of Table I by the 
method described in the earlier papers are given in Tables II and III. The 


_ reduced characteristic temperatures used for LiF and NaF have been taken 


as 330 and 175 respectively. Blackman (1935) referring to Goldmann’s 
unpublished data on the specific heats of LiF, points out that the charac- 


teristic temperature rises sharply to 690° at the low temperatures and shows 
Aa 
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TABLE I 


Frequencies of Vibration with their Griineisen Constants 























| | LiF NaF 
Frequency | 
designation Pigs Rig Me | 
Calculated v Griineisen Calculated y Griineisen 
cm.-? constant cm. constant 
ie — 
vy oo s | 307 | 2°8 246 | 2-54 
Va | +44 | 447 1-0 356 | 1-03 
(fluorine) | 
Vs | a 680 0-86 324 1-03 
(alkali) | | | 
V4 ’ 197 2-25 | 195 2-29 
(fluorine) 
Vs e 302 2-33 177 2-29 
(alkali) | 
Ve + oo | 629 0-6 H 332 0-4) 
Vz ue 351 2-2 H 280 2-53 
Vg se 382 2-0 241 2-58 
Vo aa 243 0-83 \ 150 ~—0-2 
Ven 230 | 1-34 | 122 1-5 
(Debye) | | 


no tendency of becoming constant near this value. At still lower tempera- 
tures it is liable to increase still further. If we assume that the constant 
value lies in the neighbourhood of 825° K., the reduced value comes to 
330° K. For NaF, the reduced characteristic temperature has been calcu- 
lated from the #, given by Shonka (1933) as a result of the studies on the 
intensities of X-ray reflection. 


Griineisen constants for the different frequencies have been calculated 
by the method outlined in the earlier paper (f). Owing to the occurrence 
of rox/2 in v¥{r), the repulsive parts of the force constants was split up 
into two parts, one for the neighbours and the other for the next 
nearest neighbours for the purpose of differentiation. 9/0 used in 
calculation were taken as 7°6 and 9-0 for LiF and NaF respectively. 
These values have been obtained earlier in this paper and have been 
preferred to Huggins’ values, which were those used in the case of the 
halides previously considered. If the Van-der-Waals forces are taken into 
account in determining the equilibrium conditions, the values of ro/p are 
likely to become smaller but, in the present case, there is an addi- 
tional uncertainty in the force constants due to the probable existence of 
fairly important non-central forces. The calculated values of the Griineisen 
constants are given in Table I, and the coefficients of thermal expansion 
in Tables II and III. The compressibilities are the same as were used in 
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determining the force constants. = 3 +3 has been taken equal to 2 
as in the earlier paper. For the Debye term, the Griineisen value of the 
constant as obtained from his well-known formula — has been used. 
The experimental values of the thermal expansion in the case of LiF have 
been taken from Adenstedt’s paper (/.c.) below the room temperature, and 
from Eucken and Dann@hl’s paper (/.c.) above it. The latter’s values have 
been calculated from the interpolation formula given by them. The experi- 
mental values of the specific heats are the unpublished work of Goldmann 
and have been quoted by Adenstedt. No experimenta! data, either of the 
specific heat, or of the thermal expansion, are known for NaF. Henglein 
(1925) has, however, determined the mean coefficients of expansion from 
— 79° to — 184°C. and from 0 to — 79°C. Her values for these are 
68 x 10-* and 98 x 10-* respectively, which compare well with the values 
of 64 x 10-* and 89 x 10-* evaluated graphically from the calculated values 


TABLE II 
Specific Heats and Thermal Expansion Coefficients of LiF 




















Temperature Calculated c, Experimental c, Calculated Experimental 
"Se cals./deg.gm.mol. | cals./deg. gm. mol. a x 10° aX 108 
| | 
50-9 +543 | -496 | sm 4-56 
72-5 | 1-55 1-458 | 14-8 | 12-12 
93-0 | 2-78 | 2-628 26 | 25-2 
117 | 4-17 4-093 | 40-7 41-4 
169-4 6-52 6-641 | 66 | 67-1 
193 7+28 7+452 | 72 76-0 
231-7 | 8-26 | 8-437 82- | 87-2 
271-5 | 9-00 | 9-145 89 | 96-0 
400 10-39 ss 106 112 
500 10-87 113 | 124 
600 | 11-20 | 120 137 
TABLE III 


Calculated Specific Heats and Thermal Expansion Coefficients of NaF 





Temperature | fy ax 108 
*“* | cals./deg. gm. mol. 
{ 

30 31 2:0 

50 1-507 10-7 

80 4-19 34 
100 5-50 46 
150 8-07 69 
200 | 9-45 82 

10-56 
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of Table III. Considering the approximations used, the agreement between 
the observed and the calculated values of the thermal expansion for both 
the salts must be regarded as satisfactory. 


In conclusion, the author wishes to express his grateful thanks to Pro- 
fessor Sir C. V. Raman, Kt., F.R.S., N.I., for the keen interest taken by him 
in this work. 

Summary 


The lattice spectrum, the specific heats, and the coefficients of thermal 
expansion of the lithium and sodium fluorides have been calculated on the 
basis of the Raman dynamics of crystal vibrations, bv the use of the methods 
developed by the author in the earlier papers. In the case of lithium fluoride, 
it has been found necessary to consider the mutual repulsions of the nearest 
fluorine ions in addition to those existing between the neighbouring ions 
of lithium and fluorine. An exponential form of the repulsive potential has 
been assumed, the constants of which have been evaluated from the equi- 
librium conditions of the lattice, the-infra red absorption frequency and the 
compressibility. | The coupling coefficients used by Kellermann have been 
re-evaluated to include the repulsions of the next nearest neighbours. In 
the case of lithium fluoride, for which adequate experimental data exist, a 
satisfactory agreement between the theoretical and the experimental values 
of the specific heats and the thermal expansion coefficients has been obtained. 
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7. Introduction 


A FUNDAMENTAL contribution to the theory of the solid state was made by 
Griineisen (1908) when he drew attention to the simple relation which exists, 
in the solids of relatively simple constitution, between their specific heat and 
thermal expansion. Basing his calculations (1912, 1926) on FEinstein’s 
original assumption of a single characteristic frequency of atomic vibration, 
Griineisen showed that the ratio between the specific heat and the thermal 
coefficient of expansion could be expressed as the product of the compressi- 
bility, the reciprocal of the atomic volume and a multiplier y, since known 
as the Griineisen constant. This constant is shown by the theory to be 
d log v 
dlog V’ 
quency v on the atomic volume V of the solid. With the advent of the Debye 
theory of specific heats, the Griineisen law was retained in its earliest form 
with the subsidiary assumption that the constant y was the same for all the 
frequencies in the spectrum. A theoretical proof of this assumption has 
not so far been given, but it is generally believed that the close agreement 
between the experimental values of the thermal expansion and those calcu> 
lated on the above basis for the metals and the alkali halides, is a sufficient 
proof of its validity. The experimental verification rests, however, on the 
values of the thermal expansion and the specific heat measured at moderate 
and high temperatures, where the contribution of the different frequencies 
to the thermal expansion can be sufficiently well expressed by a suitably 
assumed average value of y for all of them. That the Griineisen formula 
is not so successful in interpreting the thermal expansion of solids at very 
low temperatures is shown by the work of C. L. Lindemann (1912) who 
found a value of the ratio of the thermal expansion and the specific heat of 
rock-salt between (20-89)° K., which is very much lower than the room 
temperature value. The Griineisen assumption is, moreover, not supported 
by the experimental facts. Dr. Nayar (1941) measured the principal lattice 
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equal to — and exhibits the dependence of the characteristic fre- 





A3a 








146 Bisheshwar Dayal 


frequency of diamond at various temperatures and showed that its tempera- 
ture dependence could not be explained by the value of + used by Griineisen 
for the interpretation of the observed thermal expansion. Dr. Saxena (1944) 
has shown that Nayar’s results can be reconciled with the observed thermal! 
expansion, if it is assumed that the frequencies of diamond fall into two 
groups with different Griineisen constants. 


The forces in the case of the alkali halide crystals being exactly known, 
a direct calculation of the value of y for each frequency of the spectrum is 
possible. This has been done in the present paper for the three alkali halides, 
viz.. NaCl, KCl and KBr, whose specific heats had been previously 
studied by the author (Dayal, 1944 a, b) on the basis of the Raman theory 
of atomic vibrations in the crystal lattice (Raman, 1943). The values of y 
are found to be different for the separate frequencies indicated by that theory. 
The cubical coefficients of thermal expansion are then evaluated and found 
to be in good agreement with the experimental results. In particular, 
Lindemann’s observation of the failure of the Griineisen formula at low 
temperatures in the case of rock-salt is quantitatively explained. 


2. Griineisen’s Law in the Modified Form 


We assume that a crystal consists of 3 N tholecules. If the number of 
the atoms in a molecule is p, there are 3Np degrees of freedom which, we 
may suppose, are divided into j discrete frequencies. The statistical weight 
of a frequency v; is represented by the symbol o;, which means that 3Npo; 
degrees of freedom correspond to this frequency. The free energy F of a 


gram-molecule is given by the formula 
_ hy; 
F =4 + 3pRZo, log(1 —e - (1) 


where ¢ is the potential energy per gram-molecule of the atoms when they 
are at rest in their equilibrium positions. The equation of state can be 
obtained from the condition that the pressure P is given by 


dF 
r=-(28), ° 


V being the molecular volume. On differentiating and multiplying both 
sides by V we get 


(P + 4) V =2y9E;, (3) 


o| eS 


og ¥; 
og V’ 
quency v; and E; is the energy of oscillation corresponding to the same. 


where y; = — is the Griineisen constant pertaining to the fre- 
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This can be obtained from the Einstein energy function and is equal to 
3pR wir ie 
is assumed so that the right-hand side of the equation (3) becomes equal to 
yE. In the Debye theory the same value is obtained by assuming that y is 
the same for all the frequencies. E, however. has to be evaluated from a 
Debye energy function. If the potential energy between two atoms 


In the Einstein theory of specific heat, only one frequency 


separated by a distance r is given by — nt 5" Griineisen (1926) showed that 


A B 
¢=— v%i3 > Vis: (4) 


where m, n, a and b are constants. If we use the Griineisen notation and 
put (37) V equal to G(V), the equation of state of a solid becomes 
PV +G(V) =2 y;0; E;. (5) 
If the external pressure is zero this on to 
G(V) = ay; a; E;. 


Griineisen developed G (V) into a series of A = V - V, (cf: Roberts, 1928) 
where V, is the molecular volume at the absolute zero of temperature and 
at zero external pressure. We have then 


2 
G(V) =G(V,.) + A G'(V) +4, G" (Vy) +... (7) 
It can be shown that 
. l 
G(V.)=0 G(Y)=% 
0 | 
m+n+3 1 r (7a) 
” sa ose 7 bd | 
and G’ (V,) = 3 KV, J 


where K, is the compressibility at absolute zero of temperature and at zero 
external pressure. The value of G” (V,) has been obtained from the law of 
force (4). Substituting the value of G(V) in (6) and neglecting powers of 
A higher than the second we have 

AT m+n+3 A° 

e[-™S98) tine 
The second term in the brackets is small and can be neglected at low tempe- 
ratures so that on differentiating we get 


1f2V\_ Koy OE; 
v(st) =* = Viz 4% oF 9) 
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a is the volume coefficient of expansion and a; 7 is the contribution of 


the frequency v; to the specific heat. Equation (9) is Griineisen’s law in the 
most general form which can be used for the determination of a in conjunc- 
tion with any theory of lattice vibrations. Using Raman’s theory 
(cf. Dayal, a-b) we get for the alkali halides of NaCl structure 

a= Vea [278 (Gf) +30 (P)), 0 
where E and D are the Einstein and the Debye functions for specific heats, 
The values of s; and @p can be obtained from the papers cited above. At 
higher temperatures the second term within the brackets of the equation (8) 
has to be used so that, if we put A =K,2;y7;¢,;E; approximately we get 


x EH E (gf) +3 P(?)| 


Vo [1 - (mre tA) ee <i 


Only the correction quantity in the brackets of the denominator depends 
on the law of force assumed in (4). 


3. The Evaluation of y; 





(11) 


Griineisen considered the Finstein theory of specific heat according to 
which there is only one frequency of vibration of a solid namely, that in which 
one of the atoms is displaced about its mean position and the others remain 
at rest in their equilibrium positions. As usual the frequency can be put 


in the form : Me 

= 4 (12) 
where » is the mass and D is the combined force constant due to the dis- 
placement of the atom relative to all its neighbours. D can be split up into 
two parts D, and D, arising out of the attractive and the repulsive forces 
respectively. y is given by 

d log v lr dD, 
y= ~ log V  6Ddr (13) 

where r is the mean distance between the two neighbouring atoms of the 
lattice. Since the force constant between any two atoms is the second 
derivative of the potential energy, we have for the law of force (4), 

r D,_(m +2) D, 

D dr D 


r 9Dr_ (n + 2) D, 
and “oo ee 


! 
| 
i (14) 
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2) D 
so that y= fet4 eS Se “| (15) 
6 D 
where the proper signs have to be put for D, and D,, and D = D, + D,. 
For ionic crystals and such others for which » > m Griineisen neglected D, 


so that he obtained a constant value of y equal to : M4 . 


Since the frequency corresponding to any mode of oscillation can always 
be put in the form (12) the above method of derivation of y is quite general. 
If D, is negligible compared to D, for all the frequencies, there will be only 
one value of y given by ms 4 and the Griineisen law will be valid. Gene- 
rally, however, the term involving D, is negligible in the numerator due to 
the factor (n + 2) in the term D,, but this is not so in the denominator so 
that there is a different y for each frequency. In the case of the ionic 
crystals D, can be interpreted as a contribution of the coulomb forces to D 
and can be positive as well as negative. Its evaluation is rendered difficult 
because of the convergence difficulties. In the paper on the specific heats 
of the alkali halides (/.c.) its value was taken from the work of Kellermann 
(1940) and of Lyddane and Herzfeld (1938). Owing to the comparatively 
small contribution of this term to the numerator of (15), we have not 


ro 

Dd 

co-efficients given by Kellermann but, as in the case of equation (14) used 
by Griineisen, we have taken 

_r WD, (m+2) D,_ 3D,, 

D dr — D 2 

m being equal to one for the coulomb forces. The repulsive term D, in 

the expressions for the various frequencies of the Raman theory involves two 

(2v(r)\ 

drt) 


evaluated = rigorously from the summation terms of the coupling 


(16) 


2 2 
quantities ea c, and C2 which are respectively equal to 2 r 


and “(2 r =r where »(r) is the repulsive potential of the two 
ions (Dayal, b). These were obtained approximately from the equilibrium 
condition of the lattice and from the observed values of the compressibilities 
or the infra-red frequencies. The small contribution of the non-coulomb 
forces between the atoms other than the neighbours was thus ignored. Since 
practically the whole of c, and c, is given by the repulsive potential of Born 
and Mayer (1932) of the form 


v(r) = be-”'? (6 and p constants) 


(17) 
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4 = from this potential. We put D, in the form 
ne r) 


we can obtain 


2 
a? a + ke 


expressions for ead frequencies. Then since from (17) 
= Yv(r) sr du(r) 4 Bsuv(r) sr :D u(r) 


where k, and k, are constants. obtainable from the 


wv? oe oes r3 p or? 
we get 


5° = SACP). +2 HCP), CP) 0 


The last term within the brackets is very smal! so that it can be neglected. 


We have then 
rm, 7 
“Se"5 3 (19) 


This is of the same form as (14) if (n + 2) is replaced by ze Born and 


Mayer (/.c.) have shown that : is nearly the same for al! the alkali halides 


and can be replaced by an average value 2-9 x 108. Huggins (1937) has 
redetermined this value and has given the average value as 3:0 x 10°. We 
have used the latter value in the calculations of y. The expressions for the 
different frequencies and the values of D, and D, have been taken from our 
earlier work (6) on the specific heats of the alkali halides. For NaCl the 
numerical values of the last quantities have been taken from Lyddane and 
Herzfeld’s paper (1938), being the values used for the evaluation of the 
specific heats (Dayal, a). y; for the different frequencies are given in Table I. 
For the Debye term where the present method of evaluating y is un- 
suitable, we have used the average value given by Griineisen (1926). 
Owing to the very small contribution of the Debye term to the thermal 
expansion at all temperatures except the lowest, the choice of this value of 
y is not very important. 


The frequencies of the coupled oscillations are of the form 


mes [A +[at(1 - sitn,) + BY (ae 1] 


where A and B are the functions of the force constants, u is the reduced mass 
and m,, m, the masses of the ions. Owing to the occurrence of the + sign, 
y becomes very high for one of the coupled oscillations and very low for the 
other. In the case of the smallest transverse frequency for KC! and KBr, 
it actually becomes negative though its numerical value is still very small. 
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For KCI this value can be directly obtained since A? (1 - ee is nearly 
m, 2 
zero. In a rigorous calculation when all the forces are considered, the 


Griineisen constant of this frequency may change its value and perhaps its 
sign. But owing to its smallness it has practically no effect on the thermal 
expansion and further, due to the + sign, any change in its value is partly 
compensated by a change in the reverse direction of the y for the higher 
coupled frequency. 
TABLE I 
Theoretical Griineisen constants for the frequencies of the alkali halides 





Substance 














4. Evaluation of the Volume Coefficient of Thermal Expansion 


The values of « calculated from the equation (11) on the basis of the 
Raman theory of crystal a are given in Tables III and IV. The 


correction term ”” . me 2 .. - v; y; 3; E; depends on the inverse power law 


for the repulsive force. The woe of — according to this law is 
nearly 2 for all the alkali halides studied here. This value is not likely to 
change very much if an exponential law of force is used. We have evaluated 
the correction for NaCl with this value for ”™” a B The quantities 
involved being of the same order, the same correction has been used for the 
other two halides. Vy have been taken from Henglein’s paper (1926). 
K, for rock-salt has been calculated from the elastic constants at absolute 
zero extrapolated by Durand (1936) from his own data. For the others, 
Slater’s extrapolated values have been used (1924). The constants used in 
the calculations have been given in Table II. 


TABLE II 
Data used in the calculation of a 
| 





Substance rolp K, cm.2/dyne x 10-12 








KCl 


NaCl | , | 
KBr | | 
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TABLE III 


True calculated coefficients of cubic expansion of rock-salt 
Experimental data of Buffington and Latimer (1926) 






































Temperature a Xx 108 a Xx 10° ajc, x 108 

°K (Calculated) (Experimental) (Theoretical) 

20 2-2 7-3 

40 20-5 8-6 

80 63 = 9-2 

85 66°5 68-4 9-2 

89 69 as 9-2 

140 92-5 93-1 9-4 

200 103 106-3 9-6 

240 108 112-2 9-7 

280 112 117-5 9-8 

323 115 ws 10 
TABLE IV 
True coefficients of cubic expansion of KCI and KBr 
KCl KBr 
ee =f 
a x 10° alc, x 10° a Xx 108 alc, x 108 

20 4-2 5-9 9-7 6-4 
40 27 7°7 43 8-2 
80 68 8-1 82 8-9 
89 75 8-7 90 9-2 
140 93 9-0 101 9-2 
194 101 9-1 106 9-3 
273 109 9-4 112 9-6 
323 112 9-6 116 9-9 




















The true coefficients of expansion of NaCl as measured by Buffington 
and Latimer (1926) have been given in Table III along with the theoretical 
values. The authors have pointed out that the uncertainty of measurement 
in the total linear expansion of their specimen is of the order of 5 x 10-*cm. 
The total error due to this cause as well as to the uncertainty as regards 
temperature is of the order of 2-5 % in the neighbourhood of 300° K. Our 
calculations, however, depend on the values of Ky and r,/p each of which is 
uncertain by about 4%. The agreement between the calculated and experi- 
mental values of Table III must therefore be regarded as quite satisfactory. 

Griineisen’s theory is usually verified by finding the ratio a/c, at different 
temperatures and showing it to be a constant. We have shown the value of 
this ratio in Tables II] and IV. It can be seen that it diminishes with the 
fall of temperature in each case, though it is roughly constant above the 
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liquid-air temperature. If mean values of both a and c, over an extended 
range in this region of temperature are used, even the small differences 
would be partially smoothed out and a fairly constant value of a/c, can be 
obtained. This was actually observed by Henglein (/.c.) and her results, 
therefore, cannot be accepted as a proof of the validity of the Griineisen 
formula. A lower value of this ratio, obtained by Lindemann (J.c.) for 
rock-salt in the range (20-80)° K., is in perfect agreement with the results of 
the present calculations. For the purpose of comparison, the mean values 
of a, calculated graphically from the results of Tables III and IV, are given 
in Table V and have been compared with the experimental data. Except 
Lindemann’s result for rock-salt in the range (20-80)° K.., all the experimental 
data have been taken from Henglein’s paper. Judging from the accuracy 
of 4% claimed by the latter, the agreement between the theoretical and the 
experimental values is excellent. 








TABLE V 
Mean coefficients of expansion of the alkali halides 
NaCl | KCI | KBr 
Temperature - l ] 
“kK ax10® | ax10° | ax 10° | a x 108 a < 108 a x 10° 
(Cal.) | (Obs.) | (Cal.) | (Obs.) (Cal.) (Obs.) 
SS ES SSS: EEE SCS ETS: SURES oS S Ree OO See mene eee ae Sone ee ee ee Ee es ee ee ee 
20- 80 31 32-4 36 e 46 be 
89-194 89 92-5 91 89-7 100 101 
194-273 107 110 105 101 109 110 
273-323 113 115 110-5 110 114 118 








In conclusion, the author wishes to record his,grateful thanks to Pro- 
fessor Sir C. V. Raman, Kt., F.R.S., N.L.. for the keen interest taken by him 
in this work. 

Summary 


The Griineisen formula for the thermal expansion of the solids has been 
generalised for the case when the Griineisen constant y is not the same for 
all the frequencies of vibration. y has been evaluated for all the frequencies 
given by the Raman theory of crystal vibrations, in the case of the three 
alkali halides, viz.. NaCl, KCl and KBr whose specific heats had been pre- 
viously calculated by the author. It is found that y is, generally, different 
for each frequency. Using these values of y, the cubical coefficient of 
thermal expansion has been calculated for all these alkali halides, and shown 
to be in excellent agreement with the experimental data. It has been found 
that the ratio of the thermal expansion to the specific heat decreases gradually 
with the fall of temperature. The decrease of about 25% in the value of 
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this ratio between (20-80)°K., from the room temperature value, observed 
by Lindemann for rock-salt has thus been explained satisfactorily. The 
constant value obtained by Henglein above the liquid-air temperature has 
been shown to arise from the fact that the ratio is roughly constant in this 
region of temperature, and the differences are smoothed out by taking an 
extended range of temperature. 
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7. Introduction 


A PREVIOUS communication to these Proceedings by the author (1944) con- 
tains a detailed report on the fluorescence and absorption spectra of thirty-two 
diamonds of widely different intensities and colours of luminescence. Numerous 
sharp electronic lines of varying intensities have been found in the visible region 
in the recorded spectra. The two most prominent of these are the lines at 
4152 and 5032 A.U., and in combination with these two electronic lines, 
the vibration spectrum of the crystal lattice also appears towards longer wave- 
lengths. The appearance of the 4152 electronic line and the associated lattice 
vibrations are characteristic of blue luminescence, while the 5032 system is 
similarly associated with the yellow luminescence. Whenever the 5032 system 
is recorded with any specimen, the 4152 system is an invariable accom- 
paniment though its strength may be greater or smaller than that of the for- 
mer system. The converse is not true, in other words the 4152 system may 


and indeed often does appear without a trace of the 5032 system being 
noticed. 


The colour and intensity of luminescence may differ not only from 
diamond to diamond, but also in different parts of the same diamond: The 
luminescence patterns which arise in this way are often seen very conspicuously 
in polished cleavage plates of diamond. They have been described in detail 
by Sir C. V. Raman (1944), and their relationship to similar marked variations 
in other properties of the diamonds over their area have also been ‘remarked 
on by him. Detailed studies of the luminescence patterns (Sunanda Bai, 1944), 
the patterns of transparency in the ultraviolet (Rendall, 1944), of the 
birefringence patterns (Raman and Rendall, 1944) and of the variations of 
X-ray reflection intensity (Ramachandran, 1944) have been published. The 
existence of parallel variations in such widely different physical properties 
indicates a common physical origin, viz., a variation in the ultimate struc- 
ture of the crystal, and it thus becomes of importance to make also a spec- 
troscopic study of the luminescence patterns. In the previous investigation 
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of the author, the fluorescent light from the diamond was condensed by 
means of a cylindrical lens on the slit of the spectrograph. This procedure 
prevents the spectra from different parts of the diamond being separately 
observed. If however, a cleavage plate is exposed and an image of it when 
excited to luminescence is thrown on the slit and moved over it in a series 
of exposures, the resulting spectra exhibit the variations in the nature and 
ntensity of the light emitted by the different parts of the diamond. The 
interest of such a study is that it enables us to understand the nature of the 
luminescence patterns more fully than is possible by simple visual observa- 
tion and also to correlate the results, with those obtained in the earlier 
spectroscopic investigation. 


2. Experimental Technique 


The diamonds employed were polished cleavage plates selected from 
Sir C. V. Raman’s personal collection, to be representative of the various 
types of luminescence pattern wihch have been noticed. They have the new 
catalogue numbers N. C. 80 (D38), N.C. 82 (D235), N. C. 108 (D188), 
N.C. 110 (D190), N. C. 120 (D200) and N. C. 122 (D202), the numbering in the 
older catalogue being given within brackets. N. C. 80 was blue-fluorescent, 
N. C. 120 and N. C. 122 were yellow-luminescent, while N. C. 110, N. C. 
108 and N. C. 82 exhibited both yellow and blue luminescence. Photo- 
graphs of the luminescence patterns of all these specimens except N. C. 122 
(202) are reproduced in the plates accompanying Sir C. V. Raman’s paper 
quoted above. 


The diamonds were mounted in metal blocks with suitable apertures 
for the entrance and emergence of light. Observation of the luminescence 
pattern was made normal to the area of the plate, the illuminating beam 
being obliquely incident. The source of light for exciting fluorescence was 
a carbon are run at 220 volts with a current of 5 amperes. The light from 
the arc after passage through a water cell and a plate of Wood’s glass which 
cuts off all visible radiation, was focussed on the diamond so as to illuminate 
the whole area uniformly. The image of the diamond was carefully focussed 
on the slit of the spectrograph by means of a high-quality Sonnar lens of 5 cm. 
focal length. By moving the lens, different parts of the image were made to 
pass over the slit. The whole area of the diamond could be scanned in this 
way, a series of fluorescence spectra being obtained for each specimen studied. 
Exposures of the order of two to four hours were sufficient to record each 
spectrum in the case of the strongly luminescent diamonds N. C. 80, N. C. 
108, N. C. 110, N. C. 120 and N. C. 122, while twelve to fifteen hours were 
necessary to obtain a spectrum of N. C. 82 with satisfactory intensity. The 
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spectra were photographed on Ilford H. P2 plates by means of a Hilger 
two-prism spectrograph. All observations were made at room temperature. 
It is hoped to take up studies at liquid air temperature later. 


3. Results 


(a) Observations with diamond N. C. 80 (D38).—N. C. 80 is a striking 
example of a diamond which exhibits geometrical patterns in blue Jumines- 
cence and as visually observed, shows no trace of yellow luminescence, The 
pattern consists mainly of four pairs of intensely fluorescent bands inclined 
to each other at about 60°, interspersed by regions of non-luminescence. 
A series of four spectra were obtained for the diamond, the scanning slit 
traversing the plate so as to cut across both pairs of bands. Numerous 
bright streaks of varying intensities running throughout the length of the 
spectrum are present in every spectrogram, the pattern of streaks changing 
with the portion of the diamond scanned. Every bright streak is seen to 
correspond to a region of luminescence and the completely dark portions 
between them to regions of non-luminescence. From Fig. 1, Plate IV, where 
one of the spectra obtained is reproduced, it will be seen that the 4152 system is 
the prominent feature of the spectrum and that the 5032 system is completely 
absent. Sunanda Bai (1944) has shown that in N. C. 80 the non-luminescent 
regions transmit ultra-violet light up to 2400 A. U., whereas absorption 
commences at 3000 A. U. for the luminescent portions. 


(b) Observations with diamond N. C. 82 (D235).—This is a predominantly 
blue-luminescent diamond which however, exhibits non-luminescent areas 
towards its two ends. The region of the plate which shows blue lumines- 
cence is also traversed diagonally by a series of faint yellow bands. Eight spectro. 
grams were obtained, one of which is reproduced in Fig. 2. Four electronic 
jines at 4152, 5032, 5359 and 5895 A. U. are present in the spectrum and it will 
be noticed that the 4152 line and its subsidiary bands are more intense than 
the other three. The absolute and relative intensities of these four lines 
are different in the, different spectrograms obtained, but in the same spectro- 
gram, although the intensity of the lines vary along their length, their relative 
intensities remain constant, except where they are crossed by bright stieaks. 


The streaks exhibit a wide variation in intensity and many of them 
extend throughout the spectrum, while a number of streaks will be noticed 
to be present only in the blue, green, yellow or red regions. The ‘ blue’ 
streaks belong to the 4152 system and appear with high intensity between 
4100 and 5300 A. U. and extend weakly to 6500 A. U. The ‘ green’ streaks 
similarly belong to the 5032 system and as is to be expected are present very 
weakly in the blue also, The ‘ yellow’ and ‘ red’ streaks are independent 
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of the 4152 and 5032 systems and appear to commence from the electronic 
lines at 5359 and 5895 and extend up to 6500. It appears likely that the 
‘yellow’ and ‘ red’ streaks form two new independent systems, which may 
be called the 5359 and 5895 systems respectively. The line at 5895 deserves 
special mention. Crookes (1909) had reported a citron line in red-lumi- 
nescent diamonds, having the wave-length of the sodium doublet and which he 
attributed to the presence of sodium in diamond. The line he observed 
was evidently the 5895 line which the present investigation shows to be a 
genuine electronic line of diamond that is probably responsible for the red 
luminescence. 


Fig. 2 illustrates the different kinds of streaks obtained for N. C. 82 
and shows that the majority of them lie in the 4152 system. The lower 
half of the spectrum corresponds to the region of the diamond which contains 
non-luminescent areas interspersed by luminescent patches. It will be 
noticed that a number of ‘ blue’, ‘ green’ and ‘ yellow’ streaks are present 
in this region. The ‘ yellow’ streaks in the upper portion of the spectrum 
testify to the existence of yellow luminescence bands visually observed to 
be present in this region. 


(c) Observations with diamond N. C. 120 (D200).—As visually observed, 
this is purely yellow-luminescent showing no trace of blue fluorescence. 
The most prominent feature of its luminescence pattern is a set of four 
parallel bands. A series of six spectrograms were taken and one of these is 
shown reproduced in Fig. 3. Only the 4152 and 5032 systems appear in the 
spectrum, the 5032 system being more intense than the other. Numerous 
streaks traverse the spectra throughout the region of fluorescence, and though 
their intensities vary from point to point in the diamond, the relative inten- 
sities of the two systems are uniformly constant, every variation in the 
intensity of the 5032 system being repeated in the 4152 system. 


(d) Observations with diamond N. C. 122 (D202).— Like N. C. 120 this 
is also yellow-luminescent; its whole area, however, is covered by a series 
of parallel! yellow bands. Eight spectrograms were obtained with the 
diamond, one of which is reproduced in Fig. 4. The pattern of streaks 
is different in the different spectra obtained, but in the same spectrogram 
the 4152 and 5032 systems show identical variations in intensity. There is, 
however, one exception of a ‘ blue’ streak where the 4152 system is more 
intense than the 5032. 


(e) Observations with diamond N. C. 110(D190).—-N. C. 110 is a triangular 
shaped plate showing the mixed type of luminescence. The luminescence 
pattern consists of numerous wavy bands of blue and yellow near the apex 
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of the triangle, the rest of the plate being fairly uniformly fluorescent. A 
series of ten spectrograms were taken, the plate being scanned from the 
apex to the base. Both the 4152 and the 5032 systems appear with com- 
parable intensities. The spectra obtained in the region of uniform lumines- 
cence show little change in the relative intensities of the two systems, despite 
the many streaks that run through them. On the other hand, for spectra 
obtained near the apex, one of which is reproduced in Fig. 5, the structure 
and intensity of the ‘ blue’ and ‘ green’ streaks are seen to be quite different. 


(f) Observations with diamond N. C. 108 (D188).—N. C. 108 is a typical 
example of a diamond exhibiting both blue and yellow luminescence and 
shows a very interesting pattern consisting of an intense blue spot surrounded 
by fainter blue and yellow hexagonal rings forming a figure similar 
to that of a spider web. A series of nine spectra were obtained, showing 
an amazing number and variety of streaks. No two of the spectra 
were found to be similar, in the intensities and the pattern of the 
‘blue’, ‘ green’, ‘ yellow’ and ‘ red’ streaks obtained. The four prominent 
electronic lines are those at 4152, 5032, 5359 and 5895, with which 
are associated the four kinds of streaks mentioned above. The 4152 and 
5032 systems have comparable intensities. In many cases the streaks extend 
throughout the spectrum with approximately constant intensity. More 
interesting however, are the streaks present only in one or other of the 
four systems and whose intensities are completely independent of each other. 
In Fig. 6 is reproduced one of the spectra obtained. 


4. Significance of the Results 


It is now well established that the subsidiary bands observed in the 4152 
and 5032 systems arise from combinations of the lattice vibrations with the 
electronic lines at 4152 and 5032 A.U. Hence, the intensities of the subsi- 
diary bands should vary with that of the principal electronic lines. The 
spectroscopic study of the luminescence pattern shows this to be clearly 
the case. A close parallelism exists between the intensities of the 4152 line 
and the bands associated with it at longer wave-lengths, the variations which 
appear in the former being also clearly seen in the latter. A similar relation- 
ship is observed to exist between the 5032 line and the subsidiary bands asso- 
ciated with it. 


The variations observed in the relative intensities of the 4152 and 5032 
systems are also noteworthy. It is noticed that in blue-luminescent diamonds 
the 4152 system is present alone, the subsidiary bands of this system itself extend- 
ing to 6500 A.U. In yellow-luminescent diamonds on the other hand, both 
the systems are present, the 5032 system being the prominent feature in the 
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spectrum. The intensities of the two systems in any one diamond always 
bear a constant ratio to each other, every variation in the 5032 system appear- 
ing in the 4152 system also. In diamonds which, as visually observed, show 
both blue and yellow luminescence, the two systems appear with comparable 
intensities. But the relative and absolute intensities of the two systems vary 
not only from diamond to diamond, but in the same diamond from one 
region to the other. Among the numerous bright and faint streaks that 
cross the spectra, it is observed that some are present in both the systems 
and some in either the 4152 or 5032 systems while in the case of the third 
kind of streaks, neither their structure nor their intensities are similar in the 
two systems. In fact, fine streaks or dark regions in the 4152 system in 
many cases correspond to broad streaks in the 5032 system. 


The facts observed are comprised in the statement that diamond may 
be either non-luminescent, blue-luminescent or yellow luminescent, and 
that in diamonds exhibiting luminescence patterns the different kinds of 
behaviour may manifest themselves simultaneously in different parts of the 
same specimen lying in close juxtaposition. The differences in behaviour 
are evidently associated with differences in the ultimate structure of 
the diamond. It is particularly interesting that yellow luminescence usually 
appears in thin layers exhibiting a stratification of intensity. Besides the 
4152 and 5032 systems, we have also to consider two other and 
apparently independent kinds of luminescence appearing in the yellow 
and red regions of the spectrum respectively. They consist of the electronic 
lines at 5359 and 5895 and the subsidiary bands associated with them at longer 
wavelengths. The exact origin of these systems and their correlation with 
the 4152 and 5032 systems is a problem that demands further investigation. 


In conclusion, the author wishes to express her respectful thanks to 
Professor Sir C. V Raman for his interest and constant encouragement in 
the work. 

5. Summary 


A spectroscopic study has been made of the local variations in intensity 
and colour of the fluorescent radiation from six cleavage plates of diamond 
exhibiting blue, yellow and mixed types of luminescence. It is found that a 
close correlation exists between the variations exhibited in the intensities of 
the electronic lines at 4152 and 5032 and of their subsidiary bands, 
and also between the two systems themselves. The 4152 system is present 
alone in  blue-luminescent diamonds, while both systems appear 
in diamonds of the yellow-luminescent and mixed types. In the former, 
the 4152 system occurs weakly and a close correspondence is found to 
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exist between the intensities of the two systems. In the latter, the 
4152 and 5032 systems appear with comparable intensities. no correlation 
between them being observed except in a few cases. There is also evidence 
of the existence of two other systems of luminescent bands, associated 
with the electronic lines at 5359 and 5895 A.U. respectively. 
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